Many-body quantum ratchet in a Bose-Einstein condensate 
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We study the dynamics of a dilute Bose-Einstein condensate confined in a toroidal trap and 
exposed to a pair of periodically flashed optical lattices. We first prove that in the noninteracting case 
this system can present a quantum symmetry which forbids the ratchet effect classically expected. 
We then show how many-body atom-atom interactions, treated within the mean-field approximation, 
can break this quantum symmetry, thus generating directed transport. 
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The ratchet effect, that is, the possibility to drive di- 
rected transport with the help of zero-mean perturba- 
tions, has recently gained renewed attention due to its 
possible relevance for biological transport, molecular mo- 
tors and the prospects of nanotechnology 0, Ql • At the 
classical level, the ratchet effect can be found in periodic 
systems due to a broken space-time symmetry |3|. The 
ratchet phenomenon has also been discussed in quantum 
systems W, including the Hamiltonian limit without dis- 
sipation pj. Experimental implementations of directed 
transport range from semiconductor hctcrostructures to 
quantum dots, Josephson junctions, and cold atoms in 
optical lattices jij . 

Quantum Hamiltonian ratchets are relevant in systems 
such as cold atoms in which the high degree of quantum 
control may allow experimental implementations near to 
the dissipationless limit. Moreover, the realization of 
Bose-Einstein condensates (BECs) of dilute gases has 
opened new opportunities for the study of dynamical sys- 
tems in the presence of many-body interactions. Indeed, 
it is possible to prepare initial states with high preci- 
sion and to tune over a wide range the many-body atom- 
atom interaction. From the viewpoint of directed trans- 
port, the study of many-body quantum system is, to our 
knowledge, at the very beginning. 

In this Letter, we investigate the quantum dynamics 
of a BEG in a pair of periodically flashed optical lattices. 
We show how the interaction between atoms in the con- 
densate, studied in the mean-field approximation, can 
break the quantum symmetry present in our model in 
the noninteracting limit, thus giving rise to the ratchet 
effect. The role of noise, the validity of the mean-field 
description and the possibility to observe experimentally 
our ratchet model are discussed as well. 

We consider N condensed atoms confined in a toroidal 
trap of radius R and cross section 7rr^ , with the condition 
r R, so that the motion is essentially one-dimensional. 
The dynamics of a dilute condensate in a pair of period- 



ically kicked optical lattices at zero temperature is de- 
scribed by the Gross-Pitaevskii nonlinear equation. 
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where is the azimuthal angle, g = ?>NaR/r^ is the 
scaled strength of the nonlinear interaction (we consider 
the repulsive case, i.e., g > 0), a is the s-wave scatter- 
ing length for elastic atom-atom collisions. The kicked 
potential V{6, (j), t) is defined as 

V{0, ct>, t) = E J^i(^)<^(^ - nT) + V2{9, mt ~nT- 0], 

Vi{6)^kcose, V2{e,(j})^kcos{e-(i)), 

(2) 

where k is the kicking strength and T the period of the 
kicks. The parameters (f) € [0, 2n] and ^ S [0, T] are used 
to break the space and time symmetries, respectively. 
Note that we set ft = 1 and that the length and the 
energy are measured in units of R and fi^/mR^, with 
TO the atomic mass. The wave function normalization 
reads J^'' d9\ip{9,t)\'^ — 1 and boundary conditions are 
periodic, ip{d + 27r, t) = ■ip{0, t). 

We first consider the noninteracting case 5 = 0. Here, 
when 7^ 0, TT and ^ 7^ 0, T/2 space-time symmetries 
are broken and there is directed transport, both in the 
classical limit and, in general, in quantum mechanics 0. 
However, if we take T = Btt and ^ = 47r, then the quan- 
tum motion, independently of the kicking strength k, is 
periodic of period 2T. 

In order to prove this periodicity, it is useful to write 
the initial wave function as i/j{9,0) = '^j^Ane'xjp{in9), 
where An = ^ J^^ 'tl){9,0) exp{—in9). After free evolu- 
tion up to time t, the wave function becomes 'ip{9,t) = 

J2n "3xp (^-i^t + in9^ . If t = 47r we have ■>jj{9, t) = 

V'(6', 0) while, if t = 2ti, we obtain ^{9, t) = il;{9 + tt, 0). 
Using these relations we can easily see that the system is 



periodic with period 12n. Indeed, 

V'(6',47r+) = exp[~iVi{e)]^{0,O), 

^p{e, 67r+) = exp[-iV2{0, (j))We + TT, 47r+) 
= eM-i[V2{0, 0) - Vi(0)]}^(0 + TT, 0), 

ip{e, 107r+) = exp[-iT4(6l)]V'(e', 67r+) 
= expHy2(^,<^))^(^? + 7r,0), 

4>{e, 127r+) = exp[-iy2(6', (/>)]V'(6' + tt, 107r+) = ?A(6l, 0), 

(3) 

where ip{9,t~^) denotes the value of the wave function 
at time t just after the kick. The momentum {p{t)) — 
—i d9ip*{9,t)-^^{9,t) also changes periodically with 
period 12tt (4 kicks). Therefore, the average momentum 
Pav = limt^ooP(i) {p{t) = J J* dt' {p{t'))) is given by 



47r(p(0))+27r(p(47r+))+47r(p(67r+))+27r(p(107r+)) 
12Tr 



(P(0)) + I C - sin(0 - cj,)) \m0)\^d9. 



(4) 



In particular, for the constant initial condition ?/'(^,0) = 
1/\/2tt, which is the ground state of a particle in the 
trap, the momentum remains zero at any later time. This 
initial condition has an important physical meaning, as 
it corresponds to the initial condition for a Bose-Einstein 
condensate. 

It is therefore interesting to study the case of a BEC 
because atom-atom interactions may break the above pe- 
riodicity, and this may cause generation of momentum. 
The numerical integration of Eq. Q confirms this ex- 
pectation: as shown in Fig. ^ at g ^ the momentum 
oscillates around a mean value clearly different from zero. 
Notice that without interactions {g = 0) the momentum 
is exactly zero, so that directed transport is induced by 
the many-body atom-atom interactions. 

In Fig. El we compare the asymptotic value Pav, ob- 
tained from long numerical integrations of the Gross- 
Pitaevskii equation (dotted line with triangles), with the 
average of {p{t)) over the first 30 kicks (p(907r), continu- 
ous line with boxes). It can be seen that this short-time 
average is sufficient to obtain a good estimate of the av- 
erage momentum pav, provided that g > 0.5 It is inter- 
esting to remark that the average momentum after the 
first kicks grows monotonously with g. Therefore, the 
ratchet current provides a method to measure the inter- 
action strength in an experiment. In the inset of Fig. [3 
we show the cumulative average p{t) . For strong enough 
interactions (17 > 0.5) the convergence to the limiting 
value Pav is rather fast as we can already see from the 
main part of Fig. |21 

Since the above complete periodicity of the single parti- 
cle system = 0) is a very fragile quantum phenomenon, 
it is important to check the visibility of the ratchet effect 
when the unavoidable noise leads to a departure from 
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FIG. 1: Momentum versus time for different values of inter- 
action strength g, at fc ~ 0.74 and (j) = —n/A: g = (dashed 
line), g = 0.5 (continuous curve), g — 1 (dotted curve). 
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FIG. 2; Momentum averaged over the first 30 kicks (solid 
line with boxes) and asymptotic momentum (dotted line with 
triangles). Inset: Cumulative average p(t) as a function of 
time for different values of g. From bottom to top g = 
0.1, 0.2, 0.4, 1.0, 1.5. Parameter values: k ^ 0.74, cj> = -7r/4. 



the ideal periodic behavior. For this purpose we consider 
fluctuations in the kicking period, modeled as random 
and memoryless variations of the period between consec- 
utive kicks, with the fluctuation amplitude at each kick 
randomly drawn from a uniform distribution in the in- 
terval e,e]. We have seen that, when the the size of 
the fluctuations is e = r/200, then the ratchet current 
generated in the noninteracting case is p(907r) — —0.007. 
This value of p is much smaller than the genuine many- 
body ratchet current already shown in Fig. [3 for values 
of g > 0.2. A similar conclusion is obtained when the 
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kicks are substituted by more realistic gaussian pulses 
of width T/IO. In this case p after 30 kicks is equal to 
p(907r) = —0.01, again too small to hide or to be confused 
with the many-body ratchet effect. 

The interaction-induced generation of a nonzero cur- 
rent can be understood as follows. We approximate, for 
small values of g, the free evolution of the BEC by a 
split-operator method as in 0: 

iPie,T) w e-'^^ie-MH'^-i)\'^e-'^^iiPie,0), (5) 

~ 1 S-^ A 

where ip{9,t + At) = e~'^'^ ^p{e,t). 
In particular, we obtain 

V'(0, 47r) w exp(-i47rg|V'(6', 0)p)V'(6', 0) (6) 



and 



^{0, 27r) « eM-^^gmO, 0)|2 + \ij{9 + ir, 0){']} 
X {cos[F(6', O)]'iIj{0 + tt,0)- sm[F{e, 0)]V'(6', 0)}, 



(7) 



where we have defined F(9,0) — iTrg[ip* (9,0)ip{6 + 
7r,0) - ip{9,0)->p*{9 + 7r,0)]. Note that, in the limit 
g ^ 0, Eqs. © and iQ become i/)(6',47r) ij{d,0) and 
■0(6', 27r) = iIj{9 + tt, 0), as expected for the noninteracting 
free evolution. Using this approximation, we compute 
the evolution of the condensate for the first two kicks, 
starting from the initial condition tp{9,0) — l/V^n. We 
obtain 

V'(^^,47r+) « -^exp[~iVim exp(-z2g), 

^(^,67r+) « ^exp{-*[\^2(0,<^) ~l^i(e)]}exp(-z35) 
X {cos[r2i(6')] + sin[r2i(6l)] exp[-i2Vi(6')]} , 

(8) 

where il.i{9) = g sm{2Vi{9)). The mechanism of the 
ratchet effect is now clear: due to atom-atom interac- 
tions, the modulus square of the wave function at time 
Gtt (before the second kick) is no longer constant in 
6. Instead we have, to first order in g, \iP{6,6tt)\'^ « 
^{1-1- (?sin[4Vi(0)]}, so that the initial constant proba- 
bility distribution is modified by a term symmetric under 
the transformation 9^—9. The current after the kick 
at time t — Qir is then given by 



(p(6^+)) 



j^^ d9Vi{9,4>)mM)? 

) sin(4fccos6') = —gksm{4))Ji{Ak), 

(9) 

where Ji is the Bessel function of the first kind of index 
1. This current is in general different from zero, provided 
that 1/2(6', (f)) is not itself symmetric under 9 ~9, that 
is, when ^ 0, tt. 

In Fig. 13 we show that it is possible to control the 
direction of transport by varying the phase (jj: the cur- 
rent can be reversed simply by changing (f) — > —0. This 
current inversion can be explained by means of the fol- 
lowing symmetry considerations. The evolution of the 



gk J^^ d9 sin( 



wave-function il!{9,t) is given by Eq. After substi- 
tuting in this equation 9 —f —9, and taking into account 
that that F(-6',0,i) = V{9,~(j),t), we obtain 



4^(M) 
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gm9,t)\^ + V{9,~^,t) 



(10) 

where ip{9,t) = il}{—9,t). Therefore, if is a so- 

lution of the Gross-Pitaevskii equation, then also ip{9,t) 
is a solution, provided that we substitute </> —(j) in 
the potential V. The momentum {p(t)) of the wave- 
function 'ip{9,t) is obviously given by {p{t)) = —{p{t)), 
where {p{t)) is the momentum of ^j{9,t). This means 
that, for every ip(9,t) whose evolution is ruled by the 
Gross-Pitaevskii equation with potential V{9,(p,t), the 
wavefunction ip{9,t) evolves with exactly opposite mo- 
mentum if — > — in y. Since we start with an even 
wavefunction, tp{9,0) = ^p{—9,0) — -0(^,0), then chang- 
ing — > — changes the sign of the momentum of the 
wavefunction at any later time. 
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FIG. 3: Momentum versus time for different values of the 
parameter at « 0.74 and g = 0.5: (j> = —tt/4 (continuous 
curve), (f} = (dashed line), tf) — n/A (dotted curve). 

When studying the dynamics of a kicked BEC, it is im- 
portant to take into account the proliferation of noncon- 
densed atoms. Actually, strong kicks may lead to thermal 
excitations out of equilibrium and destroy the conden- 
sate, rendering the description by the Gross-Pitaevskii 
equation meaningless |9|, [lOj . In the following, we show 
that, for the parameter values considered in this paper, 
the number of noncondensed particles is negligible com- 
pared to the number of condensed ones, thus demonstrat- 
ing that our theoretical and numerical results based on 
the Gross-Pitaevskii equation are reliable. 

Following the approach developed in [ij (see also [loj). 
we compute the mean number of noncondensed particles 
at zero temperature as SN{t) = X^jli /o^^ ^6'|wj(6', t)p. 
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where the evohition of Vj{9,t) is determined by 



d 



H2i0,t) 

-Ht{e,t) 



(11) 

Here H^{e,t) = H{e,t) - fi{t) + gQ{t)\i^{e,t)\''Q{t), 
H{9,t) = + g\^{0,t)\'^ + V{e,(t),t) is the mean- 

field Hamiltonian that governs the Gross-Pitaevskii equa- 
tion 10), is the chemical potential {H{6,t)tlj{9,t) = 
fj.{t)'il}{0,t)), Q{t) = t-\il^{t)){ip{t)\ projects orthogonally 
to \il^(t)), and H2{e,t) = gQ{t)i;^i0,t)Q*{t). 

We integrate in parallel Eqs. ^ and The initial 

condition of the noncondensed part is obtained by diag- 
onalizing the linear operator in 11111 [lol llll | . We obtain 



u,i0,o)\ _ 1 / e + i/e 

1/4 



(12) 



where f = { p /2+2g\l(0fi)\ ) ' ^^^^ "^^^^^^ " ^Z^^'" 
initial condition of the BEG. The numerical evolution is 
performed using the split-operator method as in Eq. 10), 
with small integration steps t <^T . 

The number 5N of noncondensed particles, depending 
on the stability or instability of the condensate, grows 
polynomially or exponentially. As shown in Fig. ^ 5N 
grows polynomially at small g and exponentially for large 
g. The transition from stability to instability takes place 
a,i g — gc ~ 1.7. At 17 > gc, thermal particles proliferate 
exponentially fast, 5N ~ exp(ri), leading to a significant 
depletion of the condensate after a time td ^ \n.{N)/r. 
On the other hand, for g < g^ the exponential growth 
rate r — and the number of noncondensed particles is 
negligible for up to long times. For instance, as shown 
in Fig.H 6N « 0.2 (10) after t = 90tt (30 kicks) at g = 
0.5 (1.5), which is much smaller than the total number 
of particles N ^ 10^ - 10^ ^E^. 

Finally, we would like to comment on the experimen- 
tal feasibility of our proposal. The torus-like potential 
confining the BEC may be realized by means of optical 
billiards The kicks may be applied using a periodi- 
cally pulsed strongly detuned laser beam with a suitably 
engineered intensity, as proposed in The feasibility 
is also supported by the latest progresses in the realiza- 
tion of BECs in optical traps such as the ^^Rb BEC in 
a quasi-one-dimensional optical box trap, with conden- 
sate length 80 /im, transverse confinement 5 fim, 
and number of particles N ~ 10^ 0|. Sequences of up 
to 25 kicks have been applied to a BEC of ^^Rb atoms 
confined in a static harmonic magnetic trap, with kicking 
strength /c ~ 1 and in the quantum antiresonance case 
for the kicked oscillator model, T = 2i: Finally, 
the interaction strength g can be tuned over a very large 
range using a Feshbach resonance 0| . 
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FIG. 4: Mean number 5N of noncondensed particles versus 
time for different values of the interaction strength g: from 
bottom to top, g = 0.5, 1.5, and 2.0. Inset: 5N vs. g after 30 
kicks. Parameter values: k ^ Q.IA, cj) = — 7r/4. 
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